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Abstract 

The Heisenberg model in 3d is studied from a dual point of view. It is shown that 
it can have vortex configurations, carrying a conserved charge (U{1) symmetry). 

Vortices condense in the disordered (demagnetized) phase. A disorder pa- 
rameter (/i) is defined, dual to the magnetization (n) , which signals condensation 
of vortices, i.e. spontaneous breaking of the dual U{1) symmetry. This study 
sheds light on the procedure known as abelian projection in non abelian gauge 
theories. 
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1 Introduction 



Order disorder duality[|I], ^ plays an increasingly important role in our under- 
standing of gauge theories, specifically of QCD||^ and of its supersymmetric 
extensions 0. 

Duality exists in systems which have spatial configurations with non trivial 
topology, carrying a conserved topological charge. It consists in the possibility 
of describing the system in terms of two sets of fields: a set of local fields which 
is convenient to describe the ordered phase at weak coupling (low temperature) , 
and a dual set which is convenient at strong coupling (high temperature). The 
non zero vacuum expectation value (vev) of a local field (order parameter) ($) 
signals order. At some value of the coupling (temperature) ($) and there 
is a transition to disorder. In the disordered phase topological excitations con- 
dense, the symmetry corresponding to conservation of the topological charge is 
spontaneously broken. In terms of the dual variables the system looks ordered, 
and the non zero vev of some operator ^, (/x) carrying non zero topological 
charge is the dual order parameter. When expressed in terms of the original 
fields /i is a highly non local operator. The dual order parameter is usually 
called a disorder parameter. 

The prototype system showing duality is the 2d Ising model. The field is a 
two valued variable (T{n) — ±1, defined on the sites of a 2d square lattice. At 
low temperature there is order and (a) =^ 0, breaking the symmetry a <-> —a. 
A transition point exists where (a) — > 0; at higher temperatures the system 
looks disordered. However the same system can be described in terms of a 
dual variable a*{n*) — ±1, again on a square lattice which roughly speaking 
associates a site to each spatial link of the original lattice with the rule that 
cr* = — 1 if the values of a at the ends of the link are equal and cr* = 1 if they 
are different. 

Topologically non trivial one dimensional (spatial) exicitations exist, which 
in terms of a look as kinks (half space with ct = — 1, half space with cr = 1) and 
are highly non local. In terms of a* a kink is a configuration with a single spin 
up. When described in terms of a* (dual variable) the system is again an Ising 
model: it is self dual. 

The only change is that the Boltzman factors, (3 and (3* , of the two descrip- 
tions are related by 

Ordered region of one description correspond to disordered region of the other 
and viceversa. In the disordered phase (ct) =0, (ct*) ^ 0. (ct*) is the disorder 
parameter. 

A similar situation occurs in the 3d XY modelg (liquid He4) where 2d 
configurations (vortices) with a conserved quantum number exist. In this case 
the dual system is a Coulomb gas[|| and there is no self duality. 

In the 4d compact C/(l) gauge theory the 3d configurations with topology 
are monopoles, which condense in the disordered phase producing confinement 
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of electric charges Q. 

Not always in these systems it is possible to explicitely write the partition 
function in terms of the dual variables. 

An alternative approach is to write the disorder operator in terms of the 
original fields |^ (not their dual). This is particularly practical in numerical 
simulations. A (non local) operator is written which creates a topological exci- 
tation, and therefore carries non zero topological charge, and its vev is measured, 
in order to detect spontaneous breaking of the dual symmetry[|[ ^j. 

In the case of QCD there exists no clear idea of what the dual description 
is, except that it is presumably a gauge theory, with a gauge group independent 
of the colour group, and an interchange of roles between electric and magnetic 
charges 0, 0. 

The topological excitations which condense should have magnetic charge. 
This supports the idea that confinement is produced by dual Meissner effect in 
a dual superconducting vacuum |pT|. 

Monopoles exist in QCD in connection with any operator $ in the adjoint 
representation: they are exposed by a procedure known as abelian projectionJT2[. 

A disorder parameter can be defined, wich detects condensation of the 
monopoles for each operator $. Numerical simulations do demonstrate that 
monopoles condense inthe vacuum, independent of the choice for <i>(l3[ 

In this paper we want to study the 3d ferromagnet Heisenberg model in a 
dual way. The motivation for that, besides its intrinsic interest, is twofold 

1) Have another example of duality and a check of the construction of the 
disorder parameter. 

2) Have some insight into the abelian projection. 

The model presents topological configurations in 2 dimensions, the well known 



instantons of the 2d 0(3) a model[15, |16[ . 

We will numerically check that the disordered (high temperature phase) of 
the model is a condensate of such excitations, by use of a non local disorder 
parameter. 

The idea will prove correct and as a check finite size scaling analysis at 
the Curie point will provide a determination of the critical indices and of the 
transition temperature, which agree with the ones produced by other methods. 

To make connection with gauge theory we shall reformulate the model in 
terms of a fiber bundle. We shall also consider a gauged version of it, which is a 
2 + 1 dimensional Georgi-Glashow model. The Heisenberg model can be viewed 
as the limit of zero gauge coupling of it. The topological charge turns out to be 
the corresponding limit of t'Hooft's magnetic charge. 

The paper is organized as follows. In sect. 2 we define the model, and for- 
mulate it in terms of a fiber bundle. In sect. 3 we define the disorder parameter, 
and present its numerical evaluation, together with the determination of the 
transition point and of critical indices. Some remarks on the construction of the 
disorder parameter are contained in sect. 4. Sect. 5 describes the gauged version 
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and discusses the limit of zero gauge coupling. Sect. 6 contains some concluding 
remarks. 



2 The model. 

The partition function is 

m = / n['^"(^)] exp(-5) (2) 

with <Kl{x) the area element on the unit sphere in colour space and 
1 ^ 

^= 9/5EE(^m"(^))' ^,n{x) = [n{x + ^)-n{x)] (3) 



X 

A generic x dependent 0(3) transformation 

n{x) U{x) n{x) (4) 

leaves Z[(3] invariant, in spite of the fact that S is not invariant, since it can 
be reabsorbed in a change of the variable ^{x) which leaves the measure d^{x) 
invariant. 

The continuum version of the model is the nonlinear 0(3) cr— modelQ 

C=\{d^af (0=2 = 1) (5) 
A gauged version of the model is in the continuum 

C=\{D^df^\G^.-G^. (6) 

= 9^ - igf ■ G^^ = df,A^ - 9^1^ + gA^ A A^ 

which is a 2 + 1 dimensional Gcorgi-Glashow model, with the length of the Higgs 
field frozen to 1. 

The Heisenberg model can be considered as the limit g ^ oi the system 

(i- 

Usually an x independent frame is used in colour space, or 3 fixed unit 
vectors [i — 1, 2, 3) with 

el'-e?-^ ^.^^,=^u (7) 

However a body fixed frame (BFF) can be used[^, with unit vectors £,i{x) 
[i = 1,2,3) again obeying 

■ ij = S^J A = 6= (8) 
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and ^3(x) — n{x). 

The frame is determined modulo an x dependent arbitrary rotation around 
n{x). 



Since ^ = 1 



d^^^{x) ^ (df, A ^^{x) (9) 



or 

= (10) 
where we have defined the covariant derivative 

Df,=df,- L3f,A ^df,-if- Lu^, (11) 

Eq.(p^ also implies that 

[D,,{u;),D,{l,)]C.{x)^0 (12) 

or, by the completeness of 

[D^{uj),D,{uj)]^f-G^,^0 (13) 

LOfj, is a pure gauge 

d^uji, - d^LO^ + Lo^ Auj^ = (14) 

The geometrical meaning of the above analysis is nothing but the very definition 
of parallel transport. Eq.(|lJ) is true, apart from singularities. A consequence 
of (14) is that any field configuration n(x) can be written as a parallel transport 
from infinity to the point x of the value uq which the field has at some point at 
infinity 

n{x)—T'^xpi J T ■ Lu ^{x) dx'^^ n() = R{x)no (15) 

V means path-ordering. The choice of the path C is irrelevant as long as the field 
Lo^ is a pure gauge, or eq.(p^) is satisfied. We will show that this is not true, 
because singularities of the matrix R in eq. (p^) , which make the connection of 
the bundle non trivial. 

We now analyse the nature of such singularities. 

The model has a Noether current 

Jj^^^ = nA df'ri (16) 

corresponding to 0(3) invariance which is conserved by virtue of the equation 
of motion 

There is, however another current which is identically conserved 

Jfj. = -^Ef^afidan A dpfi (18) 
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is parallel to n, since both dafi and dfju are orthogonal to it. 

= (19) 

Also the colour invariant current 

Jf,=n - Jf^ = -^Ef^afin ■ dafi A d^n (20) 

OTT 

is identically conserved 

J,, = (21) 

The conservation law ( pT| ) is a consequence of the invariance with respect the x 
dependent rotations around n{x). The action (^ can indeed be written as 

Any rotation around n corresponds to an uj parallel to n, and leaves S invariant. 
In defining the BFF this invariance reflects in the arbitrariness by a rotation in 
the choice of ^i, ^2, being ^3 parallel to n. 

The conserved charge corresponding to ( |2l| ) is 

Q = T- / " ■ {din A d2n)(fx (22) 



1 

An 



n ■ {uji A uj2) d X 



Q is nothing but the Chern number of the 2d 0(3) cr— model, which assumes 
positive or negative integer values, 

Q = ±n (23) 
On the other hand, by use of eg. (p^) 

Q = — / fi ■ (uji A LJ2) d^x = — / {diijJ2 ■ n — 821^1 ■ n)d'^x =^ 
An J An J 

= [ (di(uj2 ■ n) - d2(uJi ■ njjd^x + 2^ [ n ■ (uJi A 0)2) d^ x 

An J An J 

and hence, since the last term equals 2Q, 

±n= f LOi ■ ftdx'' (24) 



showing that the field ciL can have a nontrivial connection. There exist singu- 
larities of ujfi where eq. (|j) is not valid. 

This can be seen by expressing |i in polar coordinates with respect to H", 
the fixed frame axes. Then one can compute uj^ getting |]l7| 

sin dd^ip 

~d^e I (25) 
— cos Od^ijj 
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Singularities can occur when cos 61 = ±1, and tp is not defined. The singularity 
is then of the form 

( ' \ 

The corresponding field is an abelian field parallel to n in the sites where n = no 

The singular field can be explicitly computed for a time independent soliton 
solution (a 2d instanton independent of xq), sitting at the origin, getting (see 
sect. 5) 

This is nothing but a Dirac string. 



3 The disorder parameter 

We will show that solitons condense in the disordered phase of the model 
{(3 < 13c), i.e. that the C/(l) symmetry eq.(|2l|) is spontaneously broken in 
the disordered phase. To do that we define a disorder parameter which is the 
vacuum expectation value (vev) of the creation operator of a soliton. 

We start by defining the creation operator of a soliton. 

Let Rq [x, y) be a singular time independent rotation which creates a soliton 
of topological charge q at site y. We will give and discuss the explicit form of 
Rq below. 

We define the lattice creation operator of a soliton at site y and time t as 
follows 

^lq{y,t) = e^Y>[~fi\^{Rq\x,y)n[x,t + l)-n{x,t)f (26) 

X 

-{n{x,t + l)-n{x,t)f I 

When computing a correlation function of /^'s the definition amounts to replace 
at the time t the time derivative term of the action [Aon(x, t)]^ (the second 
term in parenthesis of eq.(p6|)) by the first term. 
We will compute the correlator 

vixo)^{^^^qio,x'')^,q{o,o)) (27) 

i.e. the propagation of a soliton sitting in the origin from time to time xq. 
At large values of |a;o|, by cluster property one expects 

V{xa)^Aexp{~M\xa\) + {nq)^ (28) 

(fiq) ^ in the thcrmodinamic limit signals spontaneous breaking of the U{1) 
symmetry (pi|). 
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Our guess is that this spontaneous breaking is related to the phase transition, 
and that (/ig) is the disorder parameter of the system, dual to the magnetization 
in the low temperature phase. 

To check that we will compute {jjLq) numerically, and in particular we will 
study its behaviour around /3c- As a byproduct we shall determine by a finite 
size scaling analysis /3c and the critical index v of the spin correlation length, 
and (3c- 

To show that /i^ actually creates a soliton let us compute 2?(xo) of eq.(p7|) 
by use of the definition (26) of jiq 



df}exp(-/3S')/x_g(0,a;o)/ig(0,0) r„ . 

-Di..) = ^ , = ^i^^^ (29) 

dl7exp(-/3^) ^^^^ 



According to the definition ( [2^ ) S + A5 is obtained from S (eq.(|^)) by the 
replacements at the two time slices i = and t = 

t = Q (Aon(f,0))^ ^ (i?g"in(f,l) -n(f,0))^ (30) 
t^XQ {AQn{x,xo)f ^ {RZln{x,XQ + 1) - fi{x,xo)y (31) 

Z[S + AS"] can be computed by the change of variables n{x, 1) — > Rqfi{x, 1). 
As observed at the beginning of sect. 2 this change leaves the measure invariant. 
The term (AqT?)^ is restored to the primitive form, but 

(A,n)2 ^ (A,i?,n)2 

(i.e. a soliton has been created at time t = 1) and 

Aon(x, 2) {R^^n{x, 2) - n(f , 1)) 

which translates the change (^^ at t = 1. 

The construction can be repeated until time xq is reached, when the con- 
struction produces again a soliton and 

(Aon(f , xo))^ {RZ\n{x, xq + 1)- Rqn{x, xq))'^ = (Aon(x, xq))'^ 

since RZq = Rq. This shows that 'Dq{xo) actually describes a soliton propagat- 
ing from t = to t = Xq. 

Instead of the quantity p = ln(/i) will be measured Since (/^)/3=o = 1 



fi = exp 



p{x)dx 



(32) 



p{P) is shown in fig.l for different sizes of the lattice . At high /3's the integral 
defining 2?(xo) is gaussian and p can be explicitly computed. The result is 

p = ~ciL + C2 (/3 > /3c) (33) 
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with ci > 0. The comparison of this analytic computation to numerical results 
is displayed in the figure. 
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Fig.l p vs p. The perturbative evaluation at high f3 is shown by the 
dotted lines for comparison. 

(/i) is an analytic function of (3 for any finite L, and therefore it cannot 
vanish in a region of /3's without vanishing identically everywhere. Only in 
the thermodynamical limit Lee- Yang singularities are produced and (/i) can be 
identically zero for /3 > f3c, as a honest disorder parameter. 

At /9 < Pc, P tends to a finite limit asV ^ oo (fig. 2), or (/i) tends to a value 
different from zero, implying that U{\) symmetry, eg. (^l|) , is spontaneously 
broken. 

Around (3c a finite size analysis can be made to determine [3c and the critical 
index v. The argument is that, for dimensional reasons 

= (34) 

with ^ the correlation length, a the lattice spacing, L the lattice size. Near the 
critical point ^ diverges with a critical index v 

i^\Pc- f3r (35) 
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Fig. 2 p at high temperature as a function of {L size of the lattice). 
As ^ goes large, the dependence of ^ on a/^ can be neglected, and 

(M)-/(|,0) = <f>(ii/'' (/?,-/?)) (3g) 
This induces on p = dln(p) / d(3 the scaling law 

-^=F(lV^(/3,-/3)) (37) 

Corresponding to the correct value of v and f3c determinations of pL~^/'^ coming 
from lattices of different size should follow the same curve if plotted versus 
L^^^iPc — /?)■ The quality of scaling is shown in fig. 3. A best fit procedure gives 

= 0.70 ± 0.02 [0.704(6)] /3c = 0.695 ± 0.003 [0.6929(1)] (38) 

For comparison the official values[|l^, |l^, |l^ are shown in square brackets. 
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Fig. 3 Finite size scaling at the optimal values of /3c and i/. 

4 On the construction of Rq{x,y) 

We want now to analyse in detail the definition and the properties of the rotation 
Rq, introduced in the previous section (eq.(p^)). In principle the model is defined 
with the boundary condition that 

nix) ^\x\^oo no 

to have a finite action. A time independent rotation Rq{x,y) which creates a 
vortex when acting on a uniform configuration in space 

n{x,t) ^no^ (0,0,1) 

is known 

Rq{x,y) = R,{0)RAf)R;\0) (39) 

with 

e^arctg^^;^ f ^ f{r), r ^ ^/W^ (40) 
{x-y)i 

and with the boundary conditions 

/(W) = /(O) = TT (41) 
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The connection cj^ produced by Rq 



^fi — d^Rq Rq ^ 



is, in polar coordinates 



sin sin / \ / cos 

= = I cos 61 sin/ a;,. = -/' -sin 6* | (42) 

1-cos// V 



The topological charge, as given by eq.(|22|), is 

/ uje ■ fide + ^ (f Qe-ndO = cos f{0))n3{y) = 1 (43) 
47r Jci 47r Jc^ 

where the line integral is taken on a small circle Ci around the point y and a 
circle C2 at \x\ = 00. The first integral gives (1 — cos /(0))n3(?7), the second 
(1 — cos /(oo))ri3(o6). we have chosen /(O) = vr f{oo) — 0. Alternatively one 
can take /(c>o) = tt and /(O) = 0, and then the role of the two circles Ci and 
C2 is interchanged. If we act on a generic configuration n{x) 

n(x) Rqn{x) — n' 

then we can compute the variation 6Q of the topological charge produced by 

Rq 

SQ = ^ [ n{x) ■ {din A San ) d^x-^ [ n{x) ■ {din A din) d^x 
Att J Att J 

= J ^ - {diRqRq^n A d2n)d^x + ^ J n ■ {din A diRqRq^ + 

+ [ n- {diRqR:Z^nAd2RqRq^n)d^x (44) 

47r J 

Assuming the form of eq.(|39|) for Rq the result is 

[^iv)^o{y) - n(oo)a;e(cx3)] (45) 

with uje defined in eq.(||). If f{oo) = tt, /(O) = a;e(0) = 0, a;e(oo) = (0,0,2) 
and 5Q — 1. The rotation Rq changes by 1 the number of vortices. In our 
numerical simulations we have first tried a fixed b.c, nt = no = (0, 0, 1). Since, 
however, the correlation length tends to be large when approaching the critical 
point, the lattice tends to be dominated by surface effects with very little bulk 
of the system, for reasonable lattices sizes. An alternative possibility is to have 
one single line parallel to the time axis {x — y) on which n ^ uq and put the 
vortex generated by Rq at the site y. This strongly reduces the boundary effect 
and makes the lattice much bigger. A third possibility is to use plain periodic 
b.c. Then in principle eq.(p5|) as it stands from the continuum, only gives the 
correct SQ if by chance n{y) = no- More than that Rq is defined only if that is 
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true. However if we change the configuration by putting n{y) — uq the action 
only changes by ©(l/L^), the sum in eq.(Q) is not affected by that change to 
0{1/L'^). Indeed for small values of correlation length the different boundary 
conditions described above give the same result within errors. 

We finally compute the connection and the gauge field G'^^(w) for the 
configuration of a vortex. We shall refer to a configuration 

"3 = ^4^^ (46) 



^ l + \w\^\ 

with w a meromorfic funtion of degree qw = Yli fEf^- Let us put uj^ = ujj^+f^^n. 
We get 

G^.(c^) = G^. + G^, (47) 

G^^ = dt,i^^ - d^Lj^ - Lj^i A Lj^ 

Since only is defined, we will determine in such a way that G^,j{ljj) — 0. 
The solution is 

h + ih^^m^^ (48) 

w oz 

The topological charge according to eq.(|3|) is then 

Q^^i hdx' = q (49) 



47r JY 



r is the contour sum of circles around the singularities of /, i.e. around the 
poles and zeros of w. 



5 The gauge version. 

The gauge version is defined by eq.(^ in the continuum. It is a 2+1 dimensional 
Georgi-Glashow model with the size of the Higgs field frozen. On the lattice it 
becomes 

5 = /3^[i?^n(a;)]2+/3'Tr{n^4 (50) 

where 

D^j_n{x) = n{x + /t) — U^{x)n{x) (51) 

and is the plaquette. Again a body fixed frame can be defined, and with it 
w^, exactly as in sect. 2. Now 

D^n = + (-.gi*^ + lo^) ^n (52) 
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while uj^ and gA^ transform both as gauge fields. Under infinitesimal gauge 
transformation both uj^ and transform as gauge fields 

gAf, 5AgAf, + dJ (53) 

The combination cj^ — gA^j^ is covariant p7[ . In the absence of singularities 
the field can be gauged away by a transformation which brings n(x) to 
fiQ, R^^{x), where R is defined by = R{x)nQ. If the transformation is 

singular, G^v{x) is not covariant but acquires a singular additive term, parallel 
to n, Gf^uiuj) 

Gf,,^Gf,,{A) + G^,{i3) (54) 

The gauge transformation R^^{x) is called an abelian projection. Going to the 
BFF makes = and d^n — 0. The U{1) invariance corresponding to the 
rotations around n becomes in this frame a C/(l) abelian gauge invariance. The 
corresponding gauge field is — di,A^, or 



By use of eq.d5 



and 



F^,^^n-G^,^-g{A^,AA^)-n (55) 



g{Af, A 1^)3 = in • {D^n A D^n) 

= n ■ G^^ - in • {D^n A D^n) (56) 



which is the 't Hooft tensor, and is gauge invariant. Calling 

= d^F^, (57) 
j; = e.o.pF'^^ (58) 

j* is identically conserved, since in the abelian projected frame F^^^, is a curl. In 
the limit g 0, gj* coincides with the current (po|). The above analysis shows 
that the Heisenberg model can be seen as the <? ^ limit of a Higgs model. The 
abelian projection is the transformation to BFF. Its singularities depend on the 
Higgs field configurations: the gauge field for the connection is present both 
in the gauged and in the simple version of the model. 



6 Concluding remarks 

A disorder parameter (/i) has been defined for the phase transition of demag- 
netization in the 3d Heisenberg model, (fi) vanishes in the magnetized phase, 
and is non zero in the disordered phase, signalling condensation of vortices in 
the vacuum. The vortices are the instantons of the 2d version of the model. 
The critical index v and the transition temperature Pc can be determined by 
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a finite size scaling analysis, and they agree with the values obtained from the 
side of ordered phase. Duality implies a non trivial topological structure of the 
model. Vortices can be viewed as gauge singularties resulting from the abelian 
projection. In fact they have a physical role on the dynamics of the system. 
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